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THEORETICAL DETERMINATION OF WATER LOADS ON PITCHING 
HULLS AND SHOCK-MOUNTED HYDRO-SKIS 
Hy Emanuel Schnitzer 


SUMMARY 


A quasi-steady theory is developed for the unsteady plane motion of 
seaplanes with high length-beam ratios and of shock-mounted hydro-skis 
impacting on a water surface while \mdergoing pitching rotation. This 
theory is based on a dynamic -camber equivalent in which a pitching flat 
plate immersed in a stream is replaced instantaneously by a stationary 
cambered airfoil for which similar fluid particle trajectories exist at 
the boundary. Since experimental hydrodynamic data were unavailable for 
verification of the proposed theory, comparisons are made with classical 
two- and three-dimensional linearized airfoil theory for steady and 
unsteady submerged motion and with a more approximate method that neg- 
lects the rotational effect on the pressure distribution. The agreement 
with the two-dimensional unsteady oscillating airfoil theory is not very 
good because of the presence of a large unsteady circulation term but 
the three-dimensional comparisons, which include, in addition, some 
oscillating airfoil data of aspect ratio 2 , indicate fair agreement. The 
effect of rotation on the longitudinal pressure distribution for upward 
pitching is seen to broaden the stagnation peak and decrease the instan- 
taneous ratio of maximum to average pressure, whereas downward pitching 
is seen to yield the opposite result. From the comparison with the more 
approximate method that neglects rotational effects on the pressure dis- 
tribution, it is deduced that the effects of rotation might be important 
for some practical narrow-hull impacts but the more approximate theory 
could be used for the trimming shock-mounted hydro-ski cases considered. 
The loads predicted by the proposed theory for the trimming shock-mounted 
hydro-ski are less than those for the fixed-trim hydro-ski. The proposed 
theory might also be useful for calculations for a low-aspect-ratio 
pitching hydrofoil. Applications of a cambered airfoil theory to the 
determination of water-pressiire distributions on hulls with pulled up 
bows are indicated. Appendixes containing exact solutions for the pres- 
sure distribution, load, and moment on a cambered airfoil immersed in a 
stream at a finite angle of attack are included. Step-by-step computa- 
tional procedures with data-sheet headings for application of the proposed 
theory to seaplane hulls with high length-beam ratios and shock-mounted 
hydro-skis impacting on a water surface and involving rotation in pitch 
are also given in an appendix. 
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INTRODUCTION 


This paper is concerned principally with the theoretical determina- 
tion of the hydrodynamic loads and motions experienced by shock-mounted 
hydro-ski-equipped aircraft, seaplanes with high length-beam ratios, and 
other relatively narrow bodies \mdergoing unsteady planing motion on a 
water surface and involving pitching rotation. Although a water-impact 
theory for the freely trimming wide float has been developed, it is not 
believed that such a theory can be extended to cover the narrow-body case. 
Also, although much information is available for the oscillating airfoil 
covering a fairly complete range of aspect ratios from 0 to most of 
these theories are only applicable for small oscillations at small angles 
of attack. They are further restricted to deeply immersed bodies, whereas 
a hydro-ski or seaplane may have only the lower surface wetted. Another 
limitation to the use of the oscillating airfoil theory for unsteady 
planing on water arises because the wing-chord is constant, whereas the 
hydro-ski or flying -boat-hull wetted length is continually changing. 
Modifications of such airfoil theory to cover the hydrodynamic case have 
been made by Glauert and Perring (ref. l)to take into account the sur- 
face of discontinuity and by Sedov (ref. 2) to include variable body 
shape and a changing wetted length. 

In order to apply Sedov's method to nonhaxmonic motions such as 
water impacts, Fourier series solutions may be made. Such solutions are 
allowable since linearizing assmptions have been incorporated in the 
derivation. Since, however, hydro-skis operate at relatively high angles 
of attack, it is possible that such linearizations may not lead to reason- 
able approximations for actual operating conditions. This is one of the 
main reasons that these methods were not utilized in the present paper. 

In order to obtain solutions to the unsteady planing problem appli- 
cable for practical angles of attack, an approximate dynamic-camber theory 
has been developed instead for the pressure distribution on pitching 
bodies during water impact. The purposes of this paper are to present 
this theory, to corroborate its accuracy insofar as possible, and to apply 
this theory and a more approximate version thereof to solving trimming 
hull and hydro-ski Impact problems. The proposed dynamic -camber theory 
is based on the premise that a pitching flat plate deeply immersed in a 
stream may be replaced instantaneously by a stationary cambered airfoil 
for which similar flow particle trajectories exist at the boundary. It 
is believed allowable to apply this type of airfoil theory, where the 
entire foil is deeply immersed, to the hydrodynamic case where only the 
lower surface may be immersed in fluid, on the basis of Herbert Wagner's 
demonstration (ref. 3) that the longitudinal pressure distribution on 
the underside of a wing is very similar to that on the underside of a 
planing plate. The classical cambered airfoil theory is extended herein 
to include high angles of attack and finite cambers, which, when inter- 
preted in terms of the pitching flat plate, are equivalent to appreciable 
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angular velocities at high trims. This equivalence comes from the quasi- 
steady adjustment of the airfoil camber at each instant to conform approx- 
imately to the instantaneous path lines of the particles moving along the 
rotating flat plate. Aspect-ratio corrections are applied to modify the 
pressure distributions for three-dimensional flow about the cambered 
airfoil. 

In this paper an approximate theory for unsteady motion with rota- 
tion in pitch is first developed and compared with other theories and 
experimental data for oscillating wings. The load and motion equations 
for pitching hulls and hydro-skis are then developed and presented along 
with solutions for a few example cases. Several appendixes which are 
utilized in connection with the development of the theory and which pre- 
sent computational procedures for making solutions of pitching-hull and 
shock-mounted hydro-ski problems are Included. 


SYMBOIfi 


A 

B 

b 

C 


cp 


aerodynamic aspect ratio 

ratio of average pressure in transverse plane to longitudinal 
center-line pressure in plane 


beam of body 


Me 


center-of -pressure coefficient relative to step, _ - f- ■ _ 


three-dimensional lift coefficient, 

L/b 


section lift coefficient. 


P tt 2 
S cU 


I SU^ 


On 




Cn 


three-dimensional moment coefficient. 


M 


or 


M 


I ScU^ I Av^b V 


section moment coefficient. 


or 


M/b 


M/b 

£■ c^^ I 

2 2 ” 


three-dimensional hydrodynamic-force coefficient. 


N 
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section hydrodynamic-force coefficient. 


W 


beam loading coefficient, — 

pgb^ 


I 


wing chord 

aerodynamic drag force, positive in aft direction 
damping exponent 

distance along ski between pivot and step or along hull between 
center of gravity and step measured parallel to keel, positive 
when measured aft (see fig. 5) 

force taken positive in upward or aft direction 

function 

damping constant, strut compression 
damping constant, strut extension 
acceleration due to gravity 
coefficient in Blasius solution 

half height of cambered foil at maximum point (see fig. 2) 
pitching moment of inertia of trimming body 

distance between keel and center of gravity measured normal to 
keel (see fig. 5 (a)) 

spring force coefficient 

reduced-frequency parameter, cod/ 2U 


aerodynamic lift force, positive in upward direction 
length of body below undisturbed water surface (see fig. 5) 
length of body below elevated water surface (see fig. 5) 
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M 

P 

P 

q 

R 

r 

R 

s 

t 

u 


u 

V 

V 

W 

w 


pitching moment, considered positive in a nose-up direction 
distance along hydro-ski "between pivot and shock-strut attachment 
fluid pressure 

resultant velocity of fluid, |fu^ + v^ 
radius of cylindrical wing (see fig. 2) 

longitudinal distance between step or trailing edge and center 
of pressure (see fig. 5) 

real part of 

wing plan-form area 

time 

free -stream velocity at infinity for stationary body, and 

equivalent forward velocity x + « zt for moving body 

tan T 

velocity of fluid in X-direction 
resultant velocity of impacting body 
velocity of fluid in Y-direction 

weight of aircraft supported by each hydro-ski or hull 
complex potential 

horizontal displacement for determination of equivalent damber 


X forward displacement of impacting body 

y^ vertical displacement for determination of equivalent camber 

z draft or downward displacement of body normal to londisturbed 

water surfaxie 


Z = X + 
Z' = X’ 
Z" = X" 


+ iY' ^ 
+ 1Y" 

y 


coordinate axes systems for conformal transformations 
(see fig. 2) 
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a inclination of stream flow to chord of airfoil (angle of attack) 

r instantaneous circulation 

7 flight-path angle 

dc 

8 angle defining effective camber, “gu radians 


i 


velocity of hull or hydro-ski normal to keel, positive downward 


9 

K 

X 


argument of Z" 

approach parameter for free-body landing, ^ ’’’ ^o) 

length of float or hydro-ski below undisturbed water surface 
divided by mean beam 

length of float or ski below elevated water surface divided by 
mean beam 


I 

P 

T 


0L 

9(A) 


longitudinal distance along 6hord or keel measured from trailing 
edge, positive forward 

mass density of fluid 

trim of ski or hull relative to undisturbed water surface, 
positive in nose-up direction 

phase angle between angle-of -attack vector and lift vector, 
positive when angle-of -attack vector is leading 

phase angle between angle-of -attack vector and moment vector, 
positive when angle -of -attack vector is leading 

aspect-ratio correction 


0) circular frequency of oscillation 


Subscripts: 

a at pivot (see fig. 5(l>) ) 

a* at connection point of shock strut to ski (see fig. 5(li)) 

e effective value 

g gravity 
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LE 

N 

n 

0 

P 

r 

s 

X 

a 

0 , 1 , 2 ^ . . . 
a 


about leading edge 
hydrodynamic normal 
point n\miber 
at contact 

aerodynamic pitching 
hydrodynamic afterbody 
about step 

extrapolated quantity 
derivative with respect to a 
successive terms in a series 

about origin (located below midchord) (see Z-plane in fig. 2) 
maximum value 


00 infinity 

A bar over a symbol means that the symbol pertains to the mass 
center of the trimming body or to the complex conjugate in appendixes A 
and B. An asterisk denotes that rotational effects have been included. 
Dots over symbols denote the derivatives with respect to time. 


DEVEIDPMMT OF THEORY FOR UNSTEADY MOTION 
WITH ROTATION IN PITCH 
Concepts of Approach 


In order to develop an approximate theory for the oblique impact 
or unsteady planing of a relatively narrow pitching body on a water- 
surface^ the following reasoning was applied. First, for the simplified 
case of oblique impact at fixed trim, it has been shown in previous 
hydrodynamic publications that this case can be replaced at each instant 
by an equivalent planing case for which a similar load distribution 
exists. The criterion for this similarity was found to be that the hori- 
zontal velocity U of the forward intersection line of the body with the 


CONFIDENTIAL 



8 


CONFIDENTIAL 


NACA RM L56E5I. 


water surface must be equal for both 
the following sketches: 



Fixed-trim impact 


cases. This is demonstrated in 



Equivalent planing 


The effects of pitching rotation which are also taken into accoiant 
in the present paper can be divided into two parts. The first part is 
the effect of the rotational velocity in modifying the effective planing 
velocity as illustrated by the following sketches: 




Trimming impact Equivalent planing ('first approximating step) 

The second part of the rotational effect is the modification of the longi- 
tudinal distribution of load or pressure caused by the influence of body 
rotation on the flow field. Although this second effect might be handled 
by classical oscillating airfoil theory, the two-dimensional theory does 
not apply for low-aspect-ratio bodies and the available low-aspect-ratio 
theory is not believed to give closed-form solutions. Therefore, a quasi- 
steady approximation was developed based on the premise that the effect 
of pitching rotation on a flat plate could be obtained from an instanta- 
neously equivalent cambered airfoil planing on a water surface as shown 
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in the following sketches; 




Trimming impact Equivalent planing (second approximating step) 


The degree of camber 2h in the right-hand sketch, of course, depends on 
the pitching velocity. 

. - - The development and testing of this dynamic -camber concept is organ- 

ized as follows. First, the camber^ form is derived as the- quasi-steady 
equivalent of the pitching wing. Then, the longitudinal pressure distri- 
bution is determined for the entire cambered form in deeply submerged flow 
p.nd for a planing form with only the bottom side wetted. The application 
to water loads on pulled -up bows is indicated in passing. In order to 
verify the pressure distributions, loads and moments are obtained for the 
cambered form undergoing steady motion and are compared with linearized 
theory- The proposed dynamic-camber equations for the pitching wing are 
presented and tested at the zero-trim end point against classical two- 
and three-dimensional linearized airfoil theory for unsteady motion and 
with low-aspect-ratio oscillating airfoil data after incorporating the 
required aspect-ratio corrections into the proposed theory. 


Derivation of Cambered Form as Quasi-Steady Equivalent 

of Pitching Wing 

A simple quasi-steady theory for approximating the flow about a 
pitching two-dimensional flat-plate wing may be derived on the assumption 
that for relatively slow rotation the boundary for the pitching wing may 
be replaced at each instant by a different fixed-cambered boundary for 
which similar particle trajectories exist along the surface. Thus, the 
path lines traced by these particles during any small time interval would 
be similar for both cases. The fixed boundary is then defined by the 
motion of the particles during this small time interval along the flat 
surface of the plate which is rotating at some instantaneous angular 
velocity d in a field of flow having a mean translational velocity U. 
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In order to determine the form of this fixed boundary, figure l(a) which 
shows a sectional view of the rotating plate is introduced. From this 
figure the incremental horizontal distance traversed by a partical moving 
along the surface in the time dt is 

dxj^ = U dt ( 1 ) 


In the same time interval the plate has rotated through the angle da. 

The approximate downwash velocities of the particle as it moves from the 
point x-|_ to the point x^^ + dxj^ are given respectively by the equations 


y^ = -Ua - x^d 
^1 ‘^1 “ “ (^1 


( 2 ) 


The incremental change in downwash velocity expressed by dy^ should be 

identical to that for a particle moving along the hypothetical fixed- 
carabered boundary (fig. l(b)) for which the doT-mwash velocity at 
points x^ and x^ + dx^ are, respectively, written 


yi 


= u 




dx. 


dyi 




y, + dy, = U ^ + d ^ 


dx. 


dx. 


(5) 


Therefore, equating the change in downwash dy^^ for both the rotating 

plate and the stationary cambered plate and malcing use of equation (l) 
leads to the eiqjression 


dyi 

dy, = U d 3 -^ = -2Ud dt 
1 dxn 


ih) 


Equation (k) can te integrated twice with respect to after again 

making use of equation (l), to yield the shape of the equivalent stationary 
cambered airfoil y^ = ^(^l) described by the equation 

. 2 

^1 = -°“l - — (5) 
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Equation (5) may be restated in nondimens ional form by dividing by the 
half chord c/2 


y^ _ x^_a 

Tj2~ ~ Tpk 



( 6 ) 


where^ for convenience^ hjz.j'2yi is designated &. A plot of this equa- 
tion is shown in figure l(b) for a = 10° and 6 = 0.1 which are believed 
to be realistic conditions for the trimming high-length-beam-ratio sea- 
plane and for the shock-mounted trimming hydro-ski. If equation (6) is 

evaluated at the points = il (the leading and trailing edges of the 

c/2 

cambered airfoil), it is found that 



(7) 


) 


so that in figure l(b) the angles Ax = ts are seen to define the sta- 
tionary cambered form. Therefore 6 = dc/2U can be thought of as an 
effective instantaneous dynamic camber of a flat plate of chord’ c rota- 
ting at angular velocity d in a stream flawing with velocity U. 


Equation (6) is seen to represent a parabolic arc at an angle of 
attack a to the stream. It will be assumed for convenience that a cir- 
cular axe will approximate equations ( 5 ) or (6) to a reasonable degree in 
this paper for application to calculations for practical hull or hydro- 
ski-equipped seaplanes. This is permitted since the maximum deviation 
between equation (6) and a circular arc drawn through the origin and end 
points of the example airfoil of figure l(b) is of the order of 0.001c. 


Longitudinal Pressure Distribution on Equivalent 
Circular-Arc Airfoil 

In order to make solutions for the tri mm ing plate it is. first neces- 
sary to obtain the potential solution for the two-dimensional flow about 
the instantaneously equivalent circular-arc airfoil at finite angle of 
attack. This derivation is presented as appendix A and is an extension 
of the exact solution for a circular-arc airfoil at 0°. angle of attack 
given in an appendix of reference 1)-. The solution presented herein is 
based on the transformation of the circle in the z"-plane (see fig. 2) 
into the off-center circle in the Z* -plane which is then converted into 
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the circular-arc airfoil in the Z-plane by a Joukowski transform. The 
derivation was required since the pressure distribution on the lower sur- 
face of the plate must be known and it is not available from published 
thin airfoil theory. In addition, since large negative pressures are not 
believed to exist in the region of the forward water line, the concept of 
an effective length or chord Cg which extends from the rear of the air- 
foil forward to the point of zero pressure immediately ahead of the stagna- 
tion point will be applied. As stated in the introduction, justification 
for utilizing the lower surface pressure from the flow about a completely 
submerged body to represent the flow about a body planing on a surface of 
discontinuity was given by Wagner. This material is presented in figure 2k 
of reference 3 which presents a comparison of the immersed and planing 
cases up to relatively high trims. The agreement appears exceptional. 


In figure 5 are presented plots of equations (A2k) and (A 26 ) of 
appendix A showing the variation, with fractional effective chord |/cg 


of fluid pressure ratio 



on the lower svirface of two-dimensional 


cambered plates. Here i is defined as the distance from the trailing 
edge measured forward along the chord and the pressure ratio is defined 
as the pressure at a point divided by the stagnation pressure due to the 
forward velocity. A wide range of trims and cambers are covered in these 
plots. Inspection of these pressxare diagrams leads to the following two 
significant points on the behavior of a trimming hull or hydro-ski. Rela- 
tive to the fixed-trim body, the effect of rotation which is a function 
of 8 can be observed as a decrease in midchord pressure for downward 
trimming (-8) which should shift the center of pressure forward, and an 
increase in midchord pressure for upward trimming (8) which should shift 
the center of pressure aft. The effect of increasing angle of attack is 
to broaden the stagnation peak and decrease the instantaneous ratio of 
maximum pressure to average pressure while a decreasing angle of attack 
yields the opposite result. This effect may be significant for hulls 
having large bottom panels which might be designed stronger to withstand 
higher average pressures during upward pitching or having small bottom 
panels which might have to be stronger to take the high local pressures in 
downward pitching. The effects of rotation on pressure distribution may 
also be derived from classical linearized oscillating wing theory of 
reference 5* 


Application of Pressure Distributions on Cambered 
Airfoils to Hulls or Itydro-skis Having 
Longitudinally Curved Bows 

As stated previously, figure 3 presents longitudinal pressure distri- 
butions which may be directly applied in accordance with dynamic -camber 
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concepts to calculation of loads and motions of pitching hulls and hydro- 
skis with rectangular, flat bottoms. These pressure distributions may 
also be used directly for straight hulls and hydro-skis having curved-up 
bows d-uring fixed-trim Impact or planing on a water surface. In order 
to take into account the effect of aspect ratio for both of these cases 
however, a method which was devised in reference 6 is used. This method 
is described in detail in a subsequent section of this paper. When the 
plots of figure 5 are applied to bows of small upward curvature, the 
entire body may be considered to be curved longitudinally, the wetted 
chord extending from the step to the bow-water-air intersection. For 
large bow curvature, the body may be broken up into two lengths, straight 
and curved. The wetted chord of the bow arc extends from the forward 
edge of the straight section to the bow-water-air Intersection. For this 
case, the overall longitudinal pressure distribution is made up of the 
press\rre distribution for the curved section (6 finite) for the bow 
region and the distribution for the straight section wherein the whole 
body is assumed to be flat (& = O). The forward part of the distribution 
for 6=0 extending over the curved-bow region is thrown away, the 
distribution for 8 finite utilized in its place, and suitable fairing 
of the pressure curve at the body inflection point is employed. The 
normal-force coefficients required in applying the aspect-ratio correction 
may be determined from planing data obtained with the model in question 
or as a very rough approximation from the normal-force coefficients for 
the flat-plate case where the equivalent flat plate might possibly be 
represented by the wetted chord from the step to the bow-water-air inter- 
section. The trim of this chord line is the trim used in obtaining theo- 
retical or experimental normal- force coefficients for the flat-plate case. 


Loads and Moments on Cambered Airfoil 

In order to verify partially the theoretical pressure distributions 
for a cambered airfoil in deeply immersed flow as given in appendix A, 
the loads and moments obtained by Integration of the pressure equations 
for small trims might be compared with classical linearized thin-airfoil 
theory. Direct integration of the proposed parametric pressiire equa- 
tions (A 2 lj-) and (A26) in appendix A in order to obtain the loads and 
moments, however, appeared to be inconvenient so the method of Blasius 
(refs. 7 oi" 8 ) was used. Appendix B presents this solution. 

The two-dimensional lift coefficient about the circular-arc airfoil 
for finite angles of attack is derived according to the method of Blasius 
in appendix B and is given as equation (BI7). This equation is repeated 
here for convenience 


= 


_ jA _ 


fcU^ 


= 2jt(sin a cos a + cos^a tan 6) 


( 8 ) 
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which, for small angles of attack and camber, reduces to the linearized 
equation for the thin airfoil (eq. 10.8, ch. II, ref. 9) 

Cj = 2jt(a + k~) (9) 


where 4^ = tan S is a measure of the camber as shown in figure 2 of 

the present paper. The upward trimming -moment coefficient about the 
leading edge for finite trims derived in appendix B as equation (B26) is 


Cm,LE 


Mus/b 

I A" 


Sin a cos a 


cos 


^5 


+ tan 8(1 + cos^) 


( 10 ) 


which for small angles reduces to 

Si,LE “ “(I 

which is the same as the linearized equation for the thin' airfoil , 

(eq. 10.8, ch. II, ref. 9) except for the minus sign in equation (ll) 
which arises from opposite positive moment conventions between this paper 
and reference 9* 


Comparison of Dynamic-Camber Equations With 

Two- and Three-Dimensional Unsteady Linearized Airfoil Theory 

It was shown previously that the angular velocity of an airfoil in 
plane motion involving pitching rotation could be expressed as an equiv- 
alent dynamic camber defined herein as 8 = oc/2U where 8 = tan"^ 4^ 

for the equivalent stationary airfoil. If in equation (9) 4^ is con- 

sidered according to this dynamic-camber concept and dc/2U is substituted 
in its place, then equation (9) can be considered as a linearized approxi- 
mation to the two-dimensional rotational case for a flat plate. Since 
comparisons are to be made in this section with linearized theory where a 
approaches zero, the proposed theory will be altered similarly for this 
comparison with the result that tan 8 approaches 8. Equation (9) can 
therefore be rewritten as follows 

Cj* = 2rt(a + 8) = 27t(a + ^ j , (l2) 
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where the asterisk denotes the added rotational effects. Rotation is 
assumed to occur about the center of the wing chord; hence, lift forces 
due to acceleration of the virtual mass are assumed to be nonexistent. 

The equation for the moment coefficient, on the other hand, should include 
a virtual mass torque which acts in the direction opposing the angular 

2** 

acceleration and has the form - (See ref. 5*) Thus the equation 

for the two-dimensional moment coefficient about the origin including the 
effects of rotation may be written 




Jtd 

2 


2- 

jcc^g 


(13) 


In this equation the d term drops out because of symmetry about the 
origin or midchord. This term appeared as the second term on the right 

side of equation (ll) as 


A simple comparison between the linearized dynamic -camber equations 
(eqs. (12) and (l5)) and the linearized theory for the oscillating airfoil 
(ref. 5) can be made for the case of a flat-plate airfoil at 0° angle of 
attack exhibiting small sinusoidal oscillations in attitude while immersed 
in a streaming fluid. For this case the instantaneous angle of attack may 
be expressed 


a = (llj.) 

where a refers to the maximum value, i is the complex imaginary 
indicator 0) is the circular frequency of oscillation, and t is 

the time. Kie instantaneous angular velocity and acceleration are there- 
fore, respectively. 


d = = ioxt„ei<^‘t = ioxx 

dt o 


a = ^ = 

dt 


Thus equation (l2) may be restated as 


■I* = 


laxxc 
2U 

iojt . icoc 


= 2«a,e-“^il t 


(15) 

(16) 


( 17 ) 
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and, if au/2U is designated k, the reduced- frequency parameter of 
reference 5; this equation becomes 


_ 2\|i + k^ gi(ajt+tan”^k) 


(18) 


The comparable equation of Theodorsen (ref. 5) is 




= J(2F - kG)^ + (k + kF + 2G)' 


i(o)t+tan-l 


2F-kG / 


(19) 


for which the F and G functions have been evaluated and plotted in 
reference 5* (These functions are not defined in the present paper.) 

The amplitude of the lift vector of equation (18) is compared with that 
of equation (19) in figure 4(a) while the phase angles by which the lift 
leads the instantaneous angle of attack are compared in figure 4(b). 

In figure 4(a), considerable disagreement is shown between the proposed 
theory and that of reference 5 for calculating the magnitude of the lift 
vectors for the infinite-aspect-ratio case. Most of this discrepancy is 
believed to arise from the amission of certain frequency-sensitive circu- 
lation terms from equation (18) which are to be foimd in equation (l9)- 
The phase angles in figure 4(b) do not seem to be affected as severely. 


The moment coefficient equation about the origin may be derived in 
a similar manner and is stated 



(20) 


whereas the comparable Theodorsen equation is 


- F - 


, 2\2 




(I- 


G-f)‘ 


i o)t+tan' 


-1 2 




kG Z ^1 


( 21 ) 


where as before, the F and G functions are those which were tabulated 
in reference 5* The amplitude of the moment vector of equation (20) is 
compared with that of equation (2l) in figure 4(c) and the phase angles 
by which the moment leads the instantaneous angle of attank are compared 
in figure 4(d). The disagreement between the moment vectors for the 
infinite-aspect-ratio case is similar in nature and magnitude with that 
of the lift vectors. The dynamic-camber system predicts no phase 


CONFIDENTIAL 



NACA RM L56E5I 


CONFIDENTIAL 


17 


difference between the moment vector and instantaneous angle of attack, 
whereas the exact theory shows a lagging phase angle. From these compar- 
isons, it is concluded that the proposed theory should be used with 
caution for pitching bodies of large aspect ratio. 

A better picture is presented for the case of low-aspect-ratio bodies 
which are the primary consideration in this paper. The dynamic -camber 
lift equation (l2) may be further modified to include the effects of 
aspect ratio by assuming that both the. term proportional to angle of 
attack and the term proportional to angular velocity of the airfoil are 
affected in a similar way by three-dimensional flow. On the basis of 
this assumption and utilizing HeMbold's airfoil aspect-ratio correction 
for zero sweep angle (refs. 10 and ll), equation (18) becomes 

4 + Prj 



yi + k 


2 i(o)t+tan“^k) 


( 22 ) 


.The_aspe.ctrratio correction _cp(A) = — t'— is believed to be fairly 

2 + y4 + A^ 

accurate for all aspect ratios greater than 1. Equation (22) is plotted 
in figure 4(a) along with data obtained with an oscillating flat-plate 
airfoil of aspect ratio 2 (ref. 12). The three-dimensional linearized 
theory of reference 15 is also plotted for comparison. The agreement 
seems fair n.nd probably indicates that the particular circulation term 
which was so affected by frequency for the infinite-span case is probably 
small for the low-aspect-ratio case. The phase-angle comparison, plotted 
in figvire 4(b), also shows in general fair agreement for the low-aspect- 
ratio case. 

As regards the pitching moment, the Helmbold aspect-ratio correction 
cannot generally be used since the shape of the longitudinal pressure 
distribution is believed to depend on aspect ratio. Therefore, for the 
moment case, aspect ratio is taken into account as follows: The two- 

dimensional ^dynamic-camber moment equation (l3) may be extended to the 
three-dimensional case by assuming that the instantaneous longitudinal 
pressure distribution for the three-dimensional case is the same as that 
for the two-dimensional case having the same instantaneous section lift 
coefficient even though the trims are different for the two cases . This 
technique was applied to the nonrotating case in reference 6. On this 
basis, the linearized three-dimensional lift coefficient for fixed atti- 
tude may be expressed as 

Cl = 2jtcp(A)a (25) 
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vhere <p(A) may be Helmbold’s eispect-ratio correction for aspect ratios 
greater than 1. The desired two-dimensional section lift coefficient is 
obtained by dividing the actual three-dimensional lift coefficient by the 
ratio of the average pressiire in a transverse plane to the pressure at 
the intersection of that plane with the longitudinal center line. Since 
reference 6 leaves this question open, a choice of Bobyleff's coefficient 
(refs. l4 and I5) has been made. Therefore, the desired section lift 
coefficient may be expressed as 




(24) 


where 


„ _ cp(A)a 

“■e B~ 


is the effective two-dimensional angle of attack 


giving the proper longitudinal center-line pressure distribution for the 
three-dimensional case of angle of attack a. The same aspect-ratio 
correction may be used for the a term as was used for the a term as 
long as the aspect ratio exceeds 1. For the hydro-ski case, smaller 
aspect ratios are encountered so that virtual mass terms might be modi- 
fied for three dimensions by a more appropriate correction developed by 
Pabst (ref. 16) from measurements obtained on oscillating bodies submerged 
in water . 


The moment coefficient about the midchord line for a trimming flat- 
plate airfoil having an aspect ratio greater than 1 is therefore obtained 
by applying the above technique to equation (13) 


Cmp 




jtc'^g 

64U^ 


9(A) 


fjta 

I 2 


rtc^g^ 

64U^y 


0(A) 


(25) 


Since the terra proportional to a does not appear in the dynamic- 
camber moment equation (eq. (l3)) because of its symmetrical effect about 
the origin or moment axis, it is omitted here also. It is seen from 
equation (25) that, for the linearized case, the aspect-ratio correction 
may be applied to the section moment equation by direct multiplication 
as was done for the lift coefficient in equation (25). For sinusoidal 
oscillations, equation (25) becomes 



(26) 
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This moment vector is plotted in figure 4(c) along with the data of 
reference 12 and the Lawrence and Gerber theory (ref. 15). The agree- 
ment also seems fair for the low-aspect-ratio case. The phase-angle 
data plotted in figure 4(d), however, predicts no phase lag for the 
dynamic-camber equations, whereas the more exact theory and the data do. 

The above comparisons appear to substantiate use of the dynamic- 
camber concept for approximate calculation of forces and moments on low- 
aspect-ratio flat plates undergoing pitching rotation and translation 
while submerged in an infinite fluid. The extension of these concepts 
in the nonlinear form to the case of a pitching body planing, on a surface 
of discontinuity is made on the basis of the previously mentioned simi- 
larity, demonstrated by Wagner, between the lower surface pressure distri- 
butions on the submerged airfoil and the planing plate. It is further 
believed that this concept may be extended to include other body shapes - 
for example, dead rise hulls and skis - after making suitable modifica- 
tions to the design plots presented in a subsequent section of this paper. 


LOAD AND MOTION EQUATIONS FOR BODIES UNDERGOING UNSTEADY 
PLANING WITH ROTATION IN PITCH 
General Considerations 


In order to determine the loads and motions of bodies vindergoing 
combined rotation and translation through a water surface, the concept 
of an equivalent or effective planing velocity established in reference I7 
is applied and extended to include rotational effects. This velocity U 
is derived with the aid of figure 5 from which it may be observed that 
the horizontal velocity of the keel-level water-intersection point rela- 
tive to the step due to rotation is defined as 

Zt f \ 

Effective forward velocity = - ■ . (27) 

o 1X1 I 

The water rise is neglected in this calculation for simplicity since its 
effect is only noticed for short bodies as shown in reference I7 but it 
xs included in computations of wetted area and aspect ratio as in refer- 
ence 18. The effective forward velocity of an Impacting body relative 
to the undisturbed water from translational motions is 

Effective forward velocity = x + ^ (28) 

The effective forward velocity of the keel-level water-intersection point 
is obtained through the summation in the horizontal direction of the 
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effective velocities derived from the rotational motion plus the effec- 
tive velocities derived from the translational motions. Therefore, if 
equations (27) and (28) are ccmhined with the equation I = z/sin t, the 
effective planing velocity U of the keel- level water-intersection point 
with the step as a reference point is obtained as 


U = x„ + 




s tan T sin^T 


i - ly 

sin T 


(29) 


where C is defined as the downward velocity of the hull or hydro-ski 
directed normal to the keel. The instantaneous velocity in combination 
with the instantaneous values of the dynamic camber 6, the trim t,. and 
the wetted length-beam ratio determine the load and load distribution on 
the impacting body including rotational effects. 

The wetted-length— -beam ratio is determined by the substitution 

into figure 6 of A defined as the ratio of the length from step to 
keel-level water Intersection to the beam. This plot is reproduced from 
reference I8 where it was constructed from the results of a planing-data 
analysis yielding the variation of the bow water rise for a rectangular 
flat plate. 


The Instantaneous two-dimensional hydrodynamic normal-force coeffi- 


cient 


* _ 






is obtained for each value of instantaneous trim 


TT^'b 

and dynamic camber 8 = • by graphical integration of the lower- 

surface pressure distributions of figure 3 - (The expression Cg in this 
figure is replaced by A^b for the hydrodynamic case.) The resulting 
plot is presented as figure 7(a). The three-dimensional instantaneous 
hydrodynamic normal-force coefficient Cjj* is Obtained through the use 

of the three-dimensional correction plotted in figure 8 which was derived 
from experimental high-speed planing data (ref. I9). For the endpoints 
of A^,^ = 0 and A^ = 00^ the theoretical solutions given in references 20 

and 18, respectively, are plotted in figure 8 so that values of lift 
coefficient between the endpoints and the experimental length-beam ratios 
can be estimated. Since the same aspect-ratio correction is applied to 
the rotational component of the force as to the translational component 
as in equation (22), Cjj* may be substituted for Cjj in the use of 

figure 8. Therefore, the Cjg* for any trim may be found on the 

= 0 curve and is the value on the appropriate Ay curve at 

the same trim. (Refer to line of long^ dashes in fig. 8. ) 
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In order to determine the center-of-press\ire coefficient about the 
step C*p from which the instantaneous pitching moment is obtained, use 

is made of the effective trim concept mentioned in the previous section. 
This effective trim Xg is obtained by substitution of the longitudinal- 
center-line normal-force -coefficient into figure 8. The trim at 

which this value Intersects the = 0 line is the effective two- 

dimensional trim. (Refer to line of short dashes in fig. 8. ) This trim 
Tg in combination with the value of 6 previously selected may then be 
substituted into figure T(^) ’to obtain the three-dimensional dynamic- 
camber center-of-pressiare coefficient 

CSp = ^ ( 30 ) 

Fn b>v Cjj 

Figure tCi*) 'was also obtained from graphical integration of the pressure 
distributions of figure 5* 

The effects of linear and angular acceleration of the virtual mass 
are neglected in 'this analysis since these terms are believed to be small 
for the quasi-steady, low-aspect-ratio case imder consideration. The 
angular acceleration term is thought to be of the same order of magnitude 
as the term for the linear acceleration of the virtual mass normal to the 
keel, which was shown to be small for narrow bodies in reference I8. 


Load and Motion Equations for Free-to-Trim Body 
Impacting on a Water Surface 

The complete load and motion time histories for a freely trimming 
body impacting through a water s-urface may be determined by means of 
step-by-step methods of calculation. Since the instantaneous hydrodynamic 
force and moment may be computed approximately for any set of instanta- 
neous conditions by means of the dynamic-camber analysis in the previous 
section entitled "General Considerations," it is necessary only to relate 
this force and moment to the dynamics of the body to obtain the incremental 
changes in accelerations, velocities, and drafts. The entire time histo- 
ries can then be synthesized from these incremental steps. 

The required load and motion equations may be. set up with the aid of 
figure 5(a)* The hydrodynamic force is assumed to be directed normal to 
the keel since viscous forces are usually small in hydrodynamic impact or 
high-speed planing. The force equation in the normal direction is therefore 

SF = Fjj + L - W cos T = - H I (31) 

o 
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If the aerodynamic lift L, which is taken normal to the keel, is assumed 
to he equal to the normal component of the weight, L - W cos t is equal 
to zero and these terms are eliminated from equation (5l)* 1!he force 
equation for the direction parallel to the keel is omitted since for 
practical trims all forces in this direction are negligible. 

The upward trimming moment about the center of gravity is 

EM = Fjj(r - E) + Mp + Mj. = It (32) 

Although the moment Mj. arising from the hydrodynamic force on the 
afterbody is unknown at this writing, it is included to remind the reader 
that it is an extremely inq)ortant term in many practical impacts and can 
only be neglected for portions of the impact where the moment contributed 
by the afterbody is negligible. If the aerodynamic lift is assumed to be 
colinear with the normal component of the weight force (aerodynamic moment 
equals zero), then Mp goes to zero and this term is eliminated from 

equation (52). 

The effective forward velocity of the keel— level-water intersection 
in space may be obtained in a manner similar to that used in deriving 
equation (29) by vectorially adding the effective forward velocity due to 
translation of the center of gravity to the effective forward velocity 
due to rotation of the hull about the center of gravity. The resulting 
equation is 

• 

U = ^ + (E - l)r (33) 

sin T 

where ^ is obtained from integration of equation (31) and t from 
integration of equation (32). 

The draft of the step Zg is obtained from the vector addition of 
the vertical displacement of the center of gravity and the vertical dis- 
placement of the step relative to the center of gravity and may be 
expressed as 


Zg = z + E(sln T - sin TqJ + J(cos t - cos Tq) (3^) 

where the subscript 0 denotes the instant of contact between the water 
and the hull. A suggested step-by-step computational procedure for 
obtaining the load and motion time histories which utilizes the above 
equations is given in appendix C. 
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Load and Motion Equations for a Trimming ^dro-Skl Mo\mted 

on a Shock Strut and Impacting on a Water Surface 

The trimming shock-mounted hydro-ski case is handled in a manner 
similar to that for the free-to-trlm body. The equations of motion are 
derived with the aid of figure 5 (^)* The summation of forces on the 
hydro-ski normal to its keel, if the mass and inertia of the ski can be 
neglected, is 

EF = Fjj - Fa - Fa« = 0 (55) 

where the forces Fa and Fa* are considered positive when they are 

acting upward on the fuselage. The summation of the forces on the fuse- 
lage normal to the hydro-ski keel is 

EF = (L - ¥)cos T + F^ + F^, = - | t = - | ta (56) 

where the assvimption is made that the aircraft fuselage does not trim 
during the impact, and where Fa is the normal force applied at the 
pivot and Fa' is the axial force applied by the shock strut which is 

assumed to be alined normal to the keel for convenience. The fixed -trim 
assTmiption has been borne out by experiments with hydro-ski equipped air- 
planes which showed little trim change of the fuselage during the immersed 
phase of any given impact. 

The shock-stnit reaction Fa» for the general shock strut with the 
hydro-ski pivot in front as in figure 5 (b) involving damping force propor- 
tional to some power of the telescoping velocity and some arbitrary form 
of springing is 


Fa' ' k[p(T(, - T)J + g|pt|'^ (57a) 

for strut compression and 

Pa' - -K[p(t - Tq)] . o'lPtl'^ (37b) 

for strut extension where K is the spring force function which, in 
general, may depend on P^t - Tq^ and in the linear spring case reduces 

to the ordinary spring constant, G and G' are the damping coeffi- 
cients, and d is the damping exponent. For the case of the hydro-ski 
pivot behind the shock strut, G is replaced by -G and -G' is 
replaced by G' . The value of P is negative for this aft location of 
the hydro-ski pivot. 
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The upward trlanniug mcment on the hydro-ski about the pivot (a) may 
he expressed . as 


£Ma = PN(r - E) + Fa'P = laT = 0 


(58) 


This sum is equal to zero since the inertia of the hydro-ski was assumed 
to he negligible. The moment on the fuselage must he equal to zero since 
the fuselage does not trim even though fuselage inertia is finite. The 
fuselage moment equation is not given here since it would serve no useful 
purpose in view of the fixed-trim qaalification. 

The effective forward velocity of the keel— level-water Intersection 
in space given previously as equation ( 29 ) becomes in terms of the motions 
about the pivot point (a) 


U = x« + 


tan T 


, (E - Z)t 
sin T 


(59) 


The draft of the step is determined by vector addition of the vertical 
displacement of the pivot and the vertical displacement of the step rela- 
tive to the pivot and is expressed as 

Zg = Za + E^sln t - sin Tq^ (Vo) 


The velocity of the point a (and the airplane) may be found from inte- 
gration of equation ( 36 ) and the angular velocity and displacement may be 
obtained by integration of equation ( 38 ). 


DISCUSSION OF SAMPLE SOLUTIONS 


In order to demonstrate the application of the dynamic-camber theory 
to unsteady planing problems, several sample solutions have been made for 
chine-immersed Impacts. These solutions have then been compared with 
further simplified theoretical solutions for which the effect of pitching 
rotation on the pressure distribution was ignored by setting 6=0 which 
for low trims is similar to the theory of reference 21. (When 6 is set 
equal to zero for a pitching flat plate impacting through or planing on a 
water surface, the instantaneous pressure distribution on this plate is 
assumed in this paper to be the same as that for a similar plate steadily 
planing on a water surface at the same conditions of trim and draft at the 
forward velocity x which is equal to the effective forward velocity U 
of the body \mdergoing unsteady motion. Thus, the theory for 6=0 is 
temed herein as the effective forward velocity theory.) Comparisons were 
also made with still further simplified fixed-trim solutions (t = O) from 
reference I 8 in order to demonstrate the effect on the hydrodynamic 
behavior of ignoring the trimming terms. 


CONFIDENTIAL 



NACA RM L56E5I 


CONFIDENTIAL 


25 


For sinoplicity, it has been assumed in the sample solutions that 
both hull forebody and ski bottoms are rectangular flat plates, that the 
aerodynamic lift is equal to, opposite, and colinear with the weight, 
and that the hull afterbody does not contact the water. A computation by 
the proposed methods for an actual impact does not have to be restricted 
to such a high degree if estimates of aerodynamic forces, afterbody water 
loads, and conversion factors for cross-sectional shapes other than the 
flat plate are available from theory or experiment. 


Free-to-Trim Hull 

Time-history solutions are presented in figures 9 (a) and 9 (b) for 
freely trimming, narrow hulls, restrained in yaw, roll, and lateral 
motions, impacting on a water surface. (See fig. 5(a).) The solutions 
for the vertical and trimming motions for a moderate beam loading » 2, 

where are plotted in figure 9(a) with the contact conditions 

' Pgb^ 

listed under run I in table I. For this case it is evident from fig- 
ure 9(a) that the dynamic-camber solution does not differ greatly from 
the solution for 6=0 insofar as vertical acceleration of the aircraft 
is concerned but does give somewhat different rebound velocities (z at 
emergence from the water) and draft histories. The trimming velocity at 
rebound is also somewhat modified. Comparison with the fixed -trim theory 
indicates, as expected, large errors in trimming motion through use of the 
fixed-trim assumption but shows for this case that vertical accelerations 
of the center of gravity can be fairly closely estimated by the fixed- 
trim theory, whereas rebound velocity and draft history show some disagree- 
ment with trimming-theory results. Some of these results are more or less 
what would have been expected from a short extrapolation of the theoret- 
ical and experimental results for a nonchine-lmmersed, trimming float 
(ref. 22) into the moderately chine-immersed region. In this reference 
also, the effect of trimming in modifying the vertical center-of -gravity 
accelerations is shown to be small. 

As the beam loading is increased so that deeper immersions with 
larger wetted-length-beam ratio result, the effect of rotation on 

the vertical motions was expected to increase for a given pitching moment 
of inertia since 5 is proportional to This is borne out in the 

solution for a higher beam loading 16, presented in figure 9(b) 

with conditions at contact listed under run II in table I. For this case, 
figure 9(15) shows that the dynamic-camber solution departs still further 
from the solution for 6=0 although the difference in the maximum 
vertical acceleration of the center of gravity is not appreciable. The 
fixed-trim solution appears to give satisfactory results for this case 
only as a very crude approximation for the vertical motions after the 
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time of maximum acceleration. As a matter of interest, it might be 
pointed out that, in figures 9(a) and 9(h), the dynamic-camber solution 
falls between the fixed-trim solution and the solution for 8=0. 

From figures 9(a) and 9(h) and the assmption that the dynamic-camber 
theory is valid at least for estimating rotational effects, it might be 
concluded that the proposed theory might be of some value in computing 
flying -boat landing characteristics for the following cases: where the 

shape of the load-ttme history is significant j where the trimming history 
is significant in determining the initial conditions of subsequent impacts 
for seaplane porpoising calculations j or possibly for unsteady, low-aspect 
ratio hydrofoil motions for both the deeply immersed case; and the case 
where the upper surface of the foil may be unwetted. The possibility 
should not be overlooked also that, for some flying -boat configurations, 
the effect of rotation on the vertical motions might be even more signif- 
icant than that indicated in figures 9i^) ^^id 9(h) since the sample cases 
chosen were selected more or less at random and are not necessarily typi- 
cal of the worst cases. Also the effects of afterbody immersion and 
aerodynamic moments on the rotational velocity may increase the rotational 
effects on the pressure distribution. (Note that, for the sample cases 
of runs I and II, the rotational velocity at contact Tq is not very 
large.) A simple criterion for evaluation of the effects of rotation is 
the magnitude of 8 which can be estimated for different hulls and impact 
or planing conditions. 


Trimming Shock-Counted Hydro-Ski 

Time-history solutions are presented in figures 10(a) and 10(b) for 
water impacts of an aircraft equipped with pivoted hydro-skis with shock- 
strut restraint. (See fig. 5(h).) In these solutions, the following 
assumptions are made: the ski mass is negligibly small; the aircraft 

fuselage itself does not trim during the impact; the shock-strut spring 
force is proportional to the telescoping displacement; the damping force 
is proportional to the square of the telescoping velocity; and the shock- 
strut axis remains more or less normal to the hydro-ski keel. 

The solution for a moderately high hydro-ski beam loading « 16) 
with the ski pivot forward of the shock strut is plotted in figure 10(a) 
with conditions at contact listed under run III in the table of initial 
conditions . The forward pivot location type of moxanting is covered by 
the more general class of mountings for which the ski trim tends to 
decrease on contact with the water. The difference between the dynamic- 
camber solution and the solution for 8=0 for the case of figure 10(a) 
is seen to be practically insignificant. It probably would not be fair 
to compare the solutions for the trimming ski with those for the fixed- 
trim shock-mounted ski since the attachment point of the shock strut to 
the ski for the trimming case and the trim selected for the fixed-trim 
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case axe believed to be critical. In order to make comparisons between 
these two types of mounting, some rational basis of comparison taking 
into consideration the effects of differences in the large ntomber of 
existing independent variables must first be evolved. As a matter of 
interest, the fixed -trim solution (same case but with shock strut locked 
in extended position) is presented to indicate the degree of load alle- 
viation obtainable for this example case with the shock-movmted trimming 
ski. 


In figure 10(b) the solution for run IV in table I which involves 
a very high beam loading (C/^ = 132 ) and a forward location of the ski 
pivot point is presented. For this case, the departure of the dynamic- 
camber solution from the solution for 8=0 is slightly greater than 
that for figure 10 (a) as would be expected from an examination of the 
component values of 8 . The heavier loading of the ski in run IV makes 
for larger which, as for the trimming hull, results in larger values 

of 8 . 

Apparently, from figures 10(a) and 10(b) for the trimming hydro-ski 
case with the forward pivot location, the theory for 8=0 may be sub- 
stituted for the dynamic-camber theory for time-history calculations at 
least up to the time of maxim-urn acceleration. Comparison in figure 10(b) 
with the fixed-trim solution (shock strut locked in the extended position) 
again shows the load reduction achieved with the trimming shock-mounted 
ski. For this case, the reduction in maximum load is considerably smaller 
than that for the lighter beam loading case of figure 10 (a), whereas the 
ratio between the fixed and trimming times to peak load is greater for 
figure 10 (b). For both cases, the rebound velocity is believed to be 
greatly reduced for the trimming mounting. 

A further evaluation of the variation of the quantities making up 8 
indicated that, for those cases where the pivot is behind the shock strut 
or more generally where the ski trim increases on contact, the effect of 
trimming on the vertical motions might be greater. An attempt was made 
to make a niimerical solution on an electronic digital computer for this 
case which was unsuccessful due to the choice of the number of significant 
figures used in the computer to describe the design charts in this paper. 
This part of the discussion is therefore included to caution the user of 
this type of computer to make use of an extra number of significant fig- 
ures in the description of the design charts for cases with aft pivot 
locations. Although the preliminary work with the electronic digital 
computer indicated that the effect of rotation on the loads and motions 
was considerably increased for this case, further attempts were not made 
to obtain this solution since it was thought at this writing to be of less 
importance than runs III and IV. It is believed that numerical solutions 
for the case of the shock-mounted ski pivoted in the rear should not be 
exceptionally difficult where the right choices of increment size and 
number of significant figures are made. 
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CONCUISIONS • 


The proposed dynamic-camber theory has . been derived primarily for 
the analytical determination of the loads and motions of seaplanes with 
high length-beam ratios, shock-mounted hydro-skis and other bodies 
impacting or planing on a water surface while undergoing pitching rota- 
tion. Comparison of this theory with oscillaiting wing theory, with some 
low-aspect-ratio oscillating wing data, with a lower order unsteady 
planing approximation neglecting the effects of rotation on the pressure 
distribution (effective forward velocity theory), and with fixed -trim 
impact theory has led to the following conclusions: 

1. The theoretical effect of rotational velocity on the longitudinal 
pressure distribution of a pitching flat plate for an increasing angle 

of attack is to broaden the stagnation peak and decrease the instantaneous 
ratio of maximum pressure , to average pressure, whereas a decreasing angle 
of attack yields the opposite results. This effect is probably signifi- 
cant for hulls having large bottom panels which would have to be designed 
stronger to withstand high average pressures during upward pitching or 
small bottom panels which would have to be strong to take the high local 
pressures in downward pitching. 

2. For the low-aspect-ratio case, the proposed theory, with the 
exception of the sensitive moment phase angle, is in fair agreement with 
Lawrence and Gerber's three-dimensional, vmsteady, linearized airfoil 
theory and with some oscillating airfoil experiments for an aspect ratio 
of 2. For the Infinite-aspect-ratio case, because of the lack of incor- 
poration of a certain frequency-sensitive circulation term into the pro- 
posed theory, it is in relatively poor agreement with Theodorsen's exact, 
two-dimensional (infinite aspect ratio), unsteady, linearized airfoil 
theory for zero angle of attack. This circulation term is apparently 
unimportant for low aspect ratios. 

3. For freely trimming, narrow seaplanes, the dynamic -camber, 
effective-forward-velocity, or fixed-trim theories give similar values 
of maximum acceleration but, where accurate time histories of transla- 
tional and pitching motions are required, especially for the heavier beam 
loadings, a theory taking into account the effects of pitching rotation 
such as the dynamic -camber theory or some equivalent thereof is probably 
desired. 

k. For trimming hydro-skis for which the trim decreases on contact 
(pivot located forward of the shock strut), both the dynamic-camber and 
effective-forward -velocity theories give similar results. Both theories 
predict load reductions for the trimming, shock-mounted hydro-ski over 
the equivalent fixed case (shock strut locked in the extended position) 
although for the higher beam loading, the load reduction was smaller for 
the cases examined, whereas the time to peak load was increased over the 
fixed trim case. 


CONFIDENTIAL 



NACA EM L 56 E 5 I 


CONFIDEIiTIAL 


29 


5 . For low-aspect-ratio pitching hydrofoils either deeply immersed 
or planing on the surface, the dynamic-camber theory is believed to offer 
some promise. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May 2k, 1956. 
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APPENDIX A 

PRESSURE DISTRIBUTION ON A CAMBERED AIRFOIL AT 
LARGE ANGLES OF ATTACK 


This appendix is concerned with the determination of the pressure 
distribution on a cambered airfoil at large angles of attack by means of 
conformal mapping procedures. This distribution is required for obtaining 
the loads and moments on such an airfoil. 

The two-dimensional incompressible flow solution for a cambered air- 
foil at 0° angle of attack is given in appendix A of reference k. Since 
the present paper is concerned with large angles of attack, a similar 
though more general analysis is required. This solution is developed 
with the aid of figure 2 as follows: 

The exact complex potential for the flow including circulation about 
a circular cylinder of radius R for all angles of attack as represented 
by a section in the Z" plane is (see ref. 8, arts. 6.22 and 'J.12) 

w = -u(z"e"^°- + 


where a is the angle of attack between the velocity vector U of the 
fluid far from the body and the negative x"-axis and p is a constant 
denoting the circulation. 

The circle in the z"-plane may be conformally transformed into a 
cambered airfoil in the Z-plane by first mapping it on the Z' -plane where 
it is converted to another circle of radius R the center of which is 
elevated in the Y' -direction a distance h. (See fig. 2.) The converting 
function required by this process is 


f 


(Al) 


Z e 


z" = Z' - ih = Z' - ^ tan 8 (A2) 

where c/k is the distance from the origin to the intercept of the circle 

on the X*-axis and tan 8 = k^. The circle on the Z' -plane is next con- 
verted to a circular-arc airfoil of length c and maximum ordinate 2h 

in the Z-plane by the Joukowski transformation 


Z = Z* 


(cA)- 

Z' 


(A3) 
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In order to determine the pressure distribution on this circular-arc 
airfoil, the velocity distribution about it must first be known. Before 
this velocity distribution can be determined, the circulation constant T 
of equation (Al) must first be determined with the aid of the Katta con- 
dition, which predicts smooth flow at finite velocity off the trailing 
edge. Thus, the complex velocity dw/dE must be finite at the point 
Z = ^ + iO in the Z-plane which means that dw/dZ' =6 at the point 


Z' = ^ + iO in the Z' -plane. In other words, a stagnation point exists 

at the intercept of the circle and the positive X'-axis in the Z' -plane. 
Actually a stagnation point also exists at the point Z' = - ^ + iO the 

intercept of the circle and the negative X'-axis. 


The complex .velocity about the circle in the Z' -plane may be obtained 
from the complex velocity in the Z"-plane by means of the following 
equation: 


dw _ dw dZ" 
dZ' ~ dZ" dZ' ■ 


(A 4 ) 


The complex velocity in the z"-plane is obtained by taking the derivative 
of equation (Al) with respect to Z" as follows: 


d^ ^ Z"2e-^ 



(A5) 


Since dZ"/dZ' = 1 , a combination of equations (Alt-) and (A5) yields the 
complex velocity in the Z' -plane which is 


dw 

dZ' 



Z"=e-W 



(A6) 


which is identical to equation (A5). In order to evaluate this expres- 
sion at the trailing edge of the airfoil, the value of Z" at the 

stagnation point Z" “ S Re"^^ must be substituted into 

equation (a 6 ) which then becomes 


dw _ _y/g-ia R^ \ _ ir 


(A?) 
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Since the velocity at a stagnation point is equal to zero, equation (A7) 
is equated to zero and yields for the circulation constant 

r = 2UR sin (a + 6) (A8) 

If the circulation constant is inserted into equation (a 6), the complex 
velocity in the Z' -plane becomes 


dZV 


:>-ia _ 




2.,2g-ia 


2iR sln(g -f &) 
Z" 


(A9) 


The complex velocity in the Z-plane is obtained by means of the equation 


dv _ dw dZ* _ dv/dZ' 
dZ dZ* dZ dZ/dZ' 


(AlO) 


The expression for dZ/dZ' is obtained by differentiating the Joukowski 
transformation (eq. (A5)) with respect to Z' which yields 


^ = 1 . (cA)^ 

dZ' 2 


(All) 


The general complex velocity of the fluid about the cambered wing in the 
Z-plane is obtained through combination of equations (A9), (AlO), and 
(All) and gives the equation 


dw 

dZ 


-U 


,-ia _ 


R^ 2iR sin(a + 8) 

,2 -ia Z'^’ 

Z e 


Z 


2 




2 


(A12) 


The velocity at the surface of the airfoil may be determined through 
substitution of the expressions for Z" and Z' on the body into 
equation (A12). These expressions in terms of the modulus and arguments 
of Z" (R and 0, respectively) and the angle 8 denoting the camber 
magnitude are 


Z" = Re^® 


and 


(A13) 


Z' = Re^® + iR sin 8 = R(e^® + i sin 8) 


(All^) 
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The complex velocity at the airfoil surface obtained through substitution 
of equations (AI5) and {Aik) and the expression = cos 6 into equa- 
tion (AI2) is written 


dw = 
dZ 


2iUe"i® 

[sin(0 - a) + sin(a + 6)^ 

(e^® + 1 sin 6; 



(1 - 2ie^®sin 6-1 


1 


(AI5) 


This complex fluid velocity may be subdivided into its components in 
the X- and Y-directions (u and v^ respectively), by the equation 


dv 

dZ 


-u + iv 


(A16) 


The square of the resultant velocity vector q is then obtained through 
multiplication of the complex velocity equation (AI6) by its conjugate 


velocity 



-u 


iv as given by the equation 


n2 = u2 + v2 = ^ ^ ^ 

dZVdZ/ dZ dZ 


(AI7) 


This multiplication is easily accomplished through use of the following 
complex conjugate relations. If and are three complex 

quantities and if 


i = W3 


(A18) 


then 


(^) - ^^2^ " ^1^% 
and therefore equation (AI7) takes the form 


(AI9) 


= K%)(«2%)(%«3) (A20) 

Thus, q^ may be obtained from equation (AI5) by an operation similar 
to equation (A20) with the following result: 


2 = dw/^ 

^ dzW> 


= 


(21Ue~^^) (-21Ue^^) _sln(9 - a) + sln(a + 5)1^ + 1 sin 5 ) (e~^^ - i sin 5 ) 

(1 - 2ie^®sin 8 - (l + 2ie"^®sin 8 - 
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or 


2 

( 1+2 sin 9 sin 6 + sin%)^ (A 21 ) 


The distribution of the pressure p about the cambered airfoil may 
be determined by means of Bernoulli's equation which for the case of 
steady motion with viscous and gravitational forces neglected may be 
written as 



P + |pq,^ = Poo + 


(A 22 ) 


where p^ and U designate the pressure and velocity far away from the 

body. The pressure far from the body is set equal to zeroj this condition 
allows equation (A22) to be written in nondimens ional form as , 



(A23) 


1 2 

The ratio of the pressure on the body to the stagnation pressure 

is then determined throiagh combination of equations (A 21 ) and (A23) and 
is written 


= 1 - 


|pu^ 


cos^a 


B 0'\ 

2 " 2 ^ 


=0<| - I) 


(1 + 2 sin 0 sin 5 + sin%) 


(A2k) 


It is noted that this pressure ratio appears as a function of the vari- 
able 9 which is a parameter in this analysis since the position along 
the airfoil can also be obtained only as a function of 0 . The X-location 
of the pressure is determined through combination of equations (AJ), (A14), 

and the well-known relation Re^® = R(cos 0 + i sin 0)j thus. 


r 

= R •< [cos 0 + i(sin 0 + sin B)j + 


cos 9 - l(sln 9 -t- sin S.) 


cos 


^ + (sin 0 + sin B)^J 


(A25) 
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and, since ^ ^ real part of this expression becomes in 

dimensionless form 


X _ cos 9 
c 4 cos 8 



cos% 


1+2 sin 0 sin S + 

sin%/ 


(A26) 


Thus, the pressure distribution about a tvo-dimensional cambered airfoil 
deeply immersed in a fluid and at any finite angle of attack is given by 
equations (A24) and (A26) where potential flow may be assumed. 


CONFIDENTIAL 



56 


CONFIDENTIAL 


NACA RM L56E5I 


APPENDIX B 

LOADS AND MOMENTS ON A CAMBEEIED AIRFOIL OF INFINITE SPAN AT 
FINITE ANGLE OF ATTACK BY BLASIUS' THEORY 


The loads and moments on a cambered airfoil might have been obtained 
through integration of the pressiire distribution equations (A24) and (A26). 
It was believed to be simpler, however, to obtain these quantities by- 
means of Blasius' theory (ref. 7, art. 7.I or ref. 8, art. 6.4) by consid- 
ering the flow far from the body. 

According to Blasius, if the square of the complex velocity can be 
expressed by the series 

(if • • 01) 

where Hq, H-|_, H 2 . . . are independent of Z, then for large values 
of Z (far away from the body) the complex force on the body may be 
■written as 

2 + ^ = (B2) 

where D is the drag and L is the lift. The two-dimensional moment 
about the origin (located below the midchord in the Z-plane of fig. 2) 

Mjj,yb may be written as: 

^ =^{=tPiH2) (B3) 

where refers to the real part of the expression. 

In order to obtain dw/dz for the cambered wing in the appropriate 
form, equations (A2), (A3) (inverted), (A5), and (A8) are first restated 
as follows: 


Z" = Z' - ih 
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dv 

dZ" 


= -U 


-ia 


...2 


ir 


r = 2UR sin(a + 6) 


Equation (A2) transforms the symmetrically located circle in the Z"-plane 
into a displaced circle in the Z' -plane offset vertically upward a dis- 
tance h. Equation (A5) is an inverted Joukowski transform which converts 
the offset circle in the Z' -plane to a circular arc of length c and maxi- 


mum height 
radical was chosen since 


2h = ^ tan 6 


in the Z-plane. The positive sign for the 


Z and Z' must hoth become large far from the 
body. Equation (A5) gives the general complex velocity of the fluid 
about the circle of radius R in the Z" -plane with circulation F where 
the resultant velocity of fluid U at <» is inclined at an angle a to 
the X"-axis . 


The expression for the complex velocity in the Z-plane can be deter- 
mined from the complex velocity in the Z"-plane by use of the jrelation 


dw _ dw dZ 
dZ dZ" dZ 


(b 4) 


In order to obtain dZ"/dZ, Z" must first be expressed as a function 
of Zj this expression is obtained through a combination of equations (A2) 
and (A3) so that 


= I " p 


- ih 


(B5) 


If the negative root had been chosen in equation (A3) instead of the 
positive one, then, as Z approaches 00^ z" approaches 0 and this 
condition means that the flow outside the wing would be transformed to 
the inside of the circle instead of the outside of the circle for which 
the flow is defined. The derivative of Z" with respect to Z is 
therefore 


dZ" _ 1 
dZ 2 




■ (^)‘ 


(b6) 


A combination of equations (A5), (b 1|-), and (b 6) results in the expression 


dw 

dZ 


2 


211 


-(i)‘ 


-U^e 

\ 


-ia _ R^^'^ 

17 «2 i 


i£ 

Z" 


(B7) 
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In order to obtain the series expression for (dw/dZ)^ required in 
equation (B1), the series for dw/dZ will be synthesized and then squared. 
The series for the first bracketed term from equation (75^) of refer- 
ence 25 may be written aa 



where approaches 0 for large values of Z far from the body , as 

required by Blasius' solution. The series for Z" as a function of Z 
from equation (B5) and equation (755) of reference 23 may be written 


Z = Z - ih 




2 2 
vhere approaches 0. The series for Z" is therefore 


Z"^ = Z^ - 21hZ 


21 2ihf^' 


.(il.m.fgl.,.. 


(BIO) 


The series expansion of dw/dZ may be synthesized from equations (B 7 ), 
(b 8 ), (B9), and (BIQ) and may be stated as 



(Bll) 
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which, after multiplication and reduction, yields the equation for the 
first few terms which are all that axe required in this paper 


^ = -Ue“^°' 
dZ 



^ rh + UR^e^®" 



(B12) 


This expression is squared to a^ain yield the first few terms required 
by the Blasius* solution for the square of the complex velocity 


( + 2irUe ^ + guphe"^®' + 2U^^- 2U^ ^^g-2i(^ + . . . 

(BI 5 ) 

Thus, the value for of equations (Sl) and (B2) is from equations (A8) 

and (BI 3 ) 

% = 2irUe"^°' = U^ci + S) (gqg a - i sin a) (B14) 

^ COS o 


since ^ ^ combination of equations (B2) and (Bll4-) yields the 

two-dimensional lift and drag forces 


^ = jtpU^c sin(a + 5) 

b cos 5 


(BI 5 ) 


R = rtpU^c sin(a + 6) 

b cos 6 


The section lift coefficient per unit span is therefore 


7 = ' v> ' = 2jt(sin a cos a + cos^a tan 5) 

I o/ ' 


|(c)U' 


(BI 6 ) 


(BI 7 ) 


The moment is obtained in a similar manner. From equation (BI 5 ) the 
value of H 2 used in equations (b 1) and (BJ) is 


Hg = -r^ - 2Urhe"^°’ - 2U%^ + 

The real part of i^ is therefore 

2 

iH 2 = sin 2a + 2Urh sin a 


-2ia 


(BI 8 ) 


(BI 9 ) 
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The two-dimensional pitching moment about the origin expressed 
positive for increasing a is therefore obtained from combining equa- 
tions (A8), (B5), and (BI 9 ) with expressions R = ; £ and h = £ tan 6 

4 cos 6 4 

to yield the equation 


M^i _ npU^c^ 

b 5 


sin a cos a(l - tan% - tan a tan &) 


(B20) 


The section moment coefficient per unit span about the origin is therefore 

= ^ = $ sin a cos a(l - tan% - tan a tan 8) (B21) 

2 ^ 


=Mii 


In order to make comparisons with the thin airfoil low-angle derivation 
in chapter II of reference 9 (see eq. 10.8 of ref. 9)f the moment must 
be taken about the lead.ing edge. An upward pitching moment (increasing 
trim) is considered positive in this report and negative in reference 9 
and thus have opposite signs. For small angles of attack, therefore, 
the two-dimensional upward trimming moment about the leading edge may be 
stated for the present derivation 

= i(c - r) = - i(| - r + I) (B22) 

where r is the distance between the center of pressure and the trailing 
edge. Since, however, the moment about the origin is 

Mp = -I^l - rj (B25) 

the moment about the leading edge may be written 

= -L I + (B24) 

which, in combination with equations (BI 5 ) and (B20), becomes 

^ sin(a + S) - sin a cos a + sin a cos a 
b 4 L cos 5 

sin^a tan ^ 
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The upward ■trimming section momen"t coefficien"t per uni't span abou"t "the 
leading edge is therefore 


_ ^LeA _ n 

D 2„2 " 2 

2 ^ 


sin g cos a + tan 6(l+ cos^a) 


COS 


% 


(B26) 
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APPENDIX C 
METHODS OF SOLUTION 


General .- The means for making numerical solutions of the equations 
for trimming bodies and shock-mounted hydro-skis are outlined in the 
following procediires. Step-by-step processes must be resorted to, because 
the equations do not lend themselves to explicit solution. Although the 
methods suggested depend on linear extrapolation, any well-known step- 
by-step procedures may be applied. 

Both example procedures assume the immersing body to be a flat 
rectangular plate. In order to make solutions for other bottom shapes, 
the theory of references 6 and 18 or planing data obtained with the 
body of interest might be utilized in judicious modifications of the 
charts of figures 6 to 8 whereas a value of the coefficient B for average 
angles of dead rise greater than zero may be obtained from references l4, 
15, or 18. The equations in this appendix Include terms representing 
aerodynamic lift and pitching moment and hydrodynamic afterbody moment. 

In order to make fairly complete solutions by the proposed methods 
including the effects of these parameters on the overall hydrodynamic 
motions and loads, it is necessary to know the variations of the aero- 
dynamic and afterbody terms with impact geometry either theoretically 
or from experimental data.^ Both solutions presented here assume that 
both aerodynamic and hydrodynamic forces are oriented essentially normal 
to the keel and neglect the force due to acceleration of the virtual 
water mass as was done in reference I8. The shock-mounted hydro-ski 
solution assumes that the mass and moment of inertia of the hydro-ski 
may be neglected, that exponential damping exists in the shock strut, 
and that the trim of the airplane is constant during impact. The water 
rise is taken into consideration only Insofar as it affects the area 
over which the hydrodynamic force acts while the rate of change of water 
rise as a force increaser (see ref. I7) is neglected. 

The numerical solutions are made with the aid of figures 5 to 8 . 

Each time-history abscissa is subdivided into equal time increments of 
duration At where the subscript n - 1 designates the values calcu- 
lated at the end of the previous time increment and n designates the 
values being calculated. For many applications it is believed to be 
advisable to select very small increments At for the first four or 
five steps and larger increments from there on. The correct increment 
size may be established by experience acquired in making several solu- 
tions for a given problem and using different increment sizes for each 
solution. The increment size may be increased until the point is reached 
where the solutions diverge from the more accurate curve obtained with 
a very small increment size. If a small-period oscillation is present in 
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the curve, too large an increment size or too small an Increment size 
for the number of significant figures used to define the design curves 
is indicated. The first point in the time-history solution is determined 
by the initial conditions of the problem. The equations listed as com- 
putational steps outline the arrangement of computation sheets. The 
values determined from the repeated application of these equations vhen 
plotted against t give the motions of the trimming hull or the fuselage 
and shock-moimted hydro-ski as the case may be. For cases where 6 is 
small, it may be taken as zero in order to simplify the calculations. 


Computational steps for free-to-trim body with flat rectangular 
bottom .- The computational steps for a free-to-trim body with flat rec- 
tangular bottom are as follows: 


( 1 ) 


■^n = "rn-l + Tn-l(^'t) + 5 


(2) Zn = Zn-1 + Cn-l(^^^)cos %_! 


(5) Zi 


^n ^n-1 


^n ^n-l 


(At) 


(*^) ^sn = Zs(n-l) + - z^.i + E(sin - sin + J (cos Tj, - cos 


(5) An = 


■*sn 


b sin T 


n 


^n ^n-1 ^n-l^-^"^) 

(T) Tn = "Tn-l + 


(8) Un = 


Cn + (E - Anb)r 


n 


sin T 


n 


(9) ^wn “ 


(10) &n radians 

2Un 

28 . 65 f^X^b . 

&n = ,, in degrees 




n 


(11) - ^(^n^ ^n) 7(3*)) 
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( 12 ) Cjj^ = (See loop of long dashes in fig. 8) 

(«) fJn = cSn 5 V 


(lA) 



- «)°°° ^ 

w7i 


or, if wing lift may he assumed equal to the weight, 

7 - !ln! 
w 

(15) T^n = >V - 0) 

where B = 0.88 (See loop of short dashes in fig. 8) 

(• 16 ) C*pn = f(Tg„, 5n) (fig. 7(b)) 

(17) Mgn " %n^\m*^cpn 

(18) T = ^3n ~ 

“ I 

or, if the resultant aerodynamic force may be assumed to pass through 
the center of gravity, and to be equal to the weight, 

V - - 4nE 

n T 


Computational steps for trimming rectangular flat-bottomed hydro.- 
ski with shock-strut restraint .- The computational steps for trimming 
rectangular flat-bottomed hydro-ski with shock-strut restraint are as 
follows: (This method is for appreciable damping only. For cases of 
negligible damping on shock-strut extension (shock struts with dirnip 
valves), use this procedure until t is a minimum, then use the fol- 
lowing zero damping procedure until z = 0 again. ) 


(1) 


T ^(AtY 

■^n = "^n-l + Tn-l(^'t) + 


(2) z„„ = z / t \ C / ,\(At)cos T T 

an a^n-1) ^a(n-l)' ' n-1 


(5) ^an = ^a(n-l) ^ K{ 


.. 2 

C / T \COS T T (At) 


Sin 
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(4) 


b sin 


(5) = £a(n-l) + 5a(n-l)<^t> 


( 6 ) u_ - 


_ U (E - V)(^n-l * ^n-l(^4)) 


Sin 


(T) \rn " ^(^n) 6) 


(8) = ^^n-1 + {radians) = T^.^CAt)) (deg) 

( 9 ) ^n) 

(10) c]Jn) dashes in fig. 8) 

(11) Fjj^ = Cjjj^ I \m^^n 

+ fL - W)cos T 

(12) ^an = - ^S/g ^ ^ 

or, if wing lift may be assimied equal to the normal component of the 
weight, then 


(15) ^en = Vi = Oj 

where B = 0.88 (See loop of short dashes in fig. 8) 

(14) Cjpn = f(Ten. &n) (fig-T(b)) 

(15) Mgjj = Fjjn^^fn^cpn 
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(16) 


'Tn = 



^ -M^n ^Nn (^n ~ '^o)] 

C 3 P 


]l/d 


where K in this case is a general spring function of P(tj^ - Tq-), or 
for the special case of linear spring force (K = constant) and turbulent 
damping (d = 2) 


^n = ± 


-4n + EFm + ^ i^n - ^ 0 ) 


nV2 


GP' 


where and Tq are expressed in radians. These equations apply only 

for strut compression where the hydro-ski pivot is in front of the shock 
strut. For extension of the strut with this arrangement the same equa- 
tions apply with the exception that G is replaced by -G' which may 
have a different value. For the case of the hydro-ski pivot behind the 
shock strut, G is replaced by -G and -G' by +G'. Thus G changes 
sign when t goes through zero. 


(17) 


■^n = 


’’^n " ’’’n-1 


(At) 


or more accurately to prevent divergence caused by small differences 
between and ■^n-l 


‘n 


- ^^n-1 V2 

2(At) 


Computational steps for trimming rectangular flat plate hydro-ski 
with shock- strut restraint (zero damping) .- ^e computational steps for 
trimming rectangular flat plate hydro-ski with shock-strut restraint 
(zero damping) are as follows: 

( 1 ) ’’’x ~ ’’^n-5 “ ^^n-2 ^^n-l 

(Subscript x denotes extrapolation.) 


(2) zan = Za(n-1) + L(n-l)(^^)cos Xn_i 


( 5 ) Zgn “ ^s(n-l) ^a(n-l)^^^^ 

Efsin Tx - sin x^.p) 


, ^a(n-l)°°^ -^n-l^^^) , 

2 


(^) An = 


sn 

b sin X. 


CONFIDENTIAL 


NACA RM L 56 E 5 I 


CONFIDENTIAL 


47 


^an “ ^a(n-l) ^a(n-l)^^^^ 

(6) U = ~ V)('^n-1 -*• 

"" sin 

( 7 ) ?w' = f(An) (fig. 6) 

C8) &n = ^('^n-1 + radians 

5n = + ^n-l(^'t)) degrees 

Un 

(9> =Sn ' ^(^x. ^n) (fig- T(a)) 

(10) (See Loop of long dashes in fig. 8) 

(11) = Cjn I Nml' V 

,1^, V . I'to * (la - ”)““ 

w ii7i 

or if ving lift may be assumed equal to the normal component of the weighty 
then. 

^ * 

>• - 

^an " ¥ 

(13) = f(% )vvn=o) 

where B = 0.88 (See loop of short dashes in fig. 8) 

(!■*> Ojpn = fCien- Sn) (fig- 7(b)) 

(15) “ lliilNmCjpn 
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(16) T] 

for the 
For the 

where ' 

(17) +1 

( 18) T] 
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- T + ^ K-1 - ^Nn 

general spring case where K"^ is the inverse spring function, 
special case of the linear spring (K = constant), 

- T + - ^Nn 

KP 

n and Tq are expressed in radians. 


_ '*'n “ "^n-l 
At 

^ ^n ~ “^n-l 
At 
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TABLE I.- INITIAL CONDITIONS OF IMPACT 


Configuration 

Trimming hull 

Run 

I 

II 

b, ft 

10 

5 


0.88 

0.88 

Ca 

1.98 

15-9 

^ 



P ft 



e' ft 

5 

5 







P 

ft/sec 



T, slug-ft2 

1,600,000. 

1,600,000 

J, ft 

14.7 

14.7 

K 1 "h /ft. 






TJ ft /spp 


— 




Vq, ft/sec 

135 

155 

W, lb 

123,500 

125,500 










0 

0 

•7 ft / - * 



^ 3 .q^ ^ 



^Q, ft/sec 

22.2 

22.2 

3 , deg 

0 

0 

7o, deg 

9.48 

9.48 

t ft. / ROP 



;ao^ ' 



t^, ft/sec 

54.9 

54.9 

to, 

0 

0 

ft/sec^ 



1 C 

0.592 

0.592 

p, slug-ft^ 

1.958 

1.958 

Tq> deg . 

5.5 

5.5 

Tq, radian/sec 

-0.087 

-0.087 

Tq, radian/sec^ 

0 

0 


Shock-mounted 
trimming ski 


III 

2 

0.88 

16.26 

2 

10 

15 

412 


IV 

1 

0.88 

132 

2 

10 

15 

412 


7,500 

223.5 


7,500 

223.5 


8,250 

0 


8,250 
0 


17 


17 


0 

6.06 

57.9 


0 

6.06 

57.9 


0 

2.23 

1.958 

15 

0 

0 


0 

2.23 
1.938 
15 
0 
0 
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.9 1.0 


Fractional effective chord, — 
c_ 


Fractional effective chord . — 

' Ce 


Figure 3.- Variation of lower surface -pressure coefficients with 
fractional effective chord for two-dimensional cambered plates 
having a wide range of camber 6. 
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Figure 5*- Concluded, 
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Figure 4 ,- Concluded 




X 


(b) Shock-moimted hydro-skis. 

Figure Geometric relations of trimming bodies during water impact. 
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Figure 7.- Variation of two-dimensional normal-force and center-of -pressure coefficients on the 
lower siorface of a cambered plate with trim t and effective camber 6. 
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6l 



10 

Trim,T, deg 


Figure 8.- Variation of planing normal -force coefficient on the level 
surface of a flat rectangular plate with trim and wetted-length- 
beam ratio \f. 
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-Proposed theory (6 varies) 

-Effective forward velocity theory (6=0) 
-Fixed trim theory (t = O) 


k 




Time t, sec 


Time t, sec 


(a) Run I. 

Figure 9*- Comparison of approximate theories for computing motion-time 
histories of freely trimming, narrow flying boats dtiring water impact 
For initial conditions see table I. 
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(a) Run III, forward pivot location. 

Figure 10.- Comparison of approximate theories for computing motion-time 
histories during water impacts of aircraft equipped with trimming, 
shock-strut-mounted hydro-skis . For initial conditions see tahle I . 
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Dynamic -camber theory (6 varies) 

c\j — Effective forward velocity theory (6 = 0) 

ffl — — — Fixed -trim theory (t = O) 

GO 



c\j — Effective forward velocity theory (6 = 0) 

ffl — — — Fixed -trim theory (t = O) 

GO 
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